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Moment closures of the Vlasov-Ampere system, whereby higher moments are represented as functions of lower 
moments with the constraint that the resulting fluid system remains Hamiltonian, are investigated by using 
water-bag theory. The link between the water-bag formalism and fluid models that involve density, fluid 
velocity, pressure and higher moments is established by introducing suitable thermodynamic variables. The 
cases of one, two and three water-bags are treated and their Hamiltonian structures are provided. In each case, 
we give the associated fluid closures and we discuss their Casimir invariants. We show how the method can be 
extended to an arbitrary number of fields, i.e., an arbitrary number of water-bags and associated moments. 
The thermodynamic interpretation of the resulting models is discussed. Finally, a general procedure to derive 
Hamiltonian V-field fluid models is proposed. 


I. INTRODUCTION 

Due to their high temperature, many plasmas, such as 
the ones encountered in the core of tokamaks or in the 
magnetosphere, can be considered as collisionless. Conse¬ 
quently, they may be well described by kinetic equations 
such as the Vlasov-Maxwell system where the particle dy¬ 
namics is described by a distribution function /(x, v, t) 
defined on a six-dimensional phase space (x, v) with x 
being the position and v the velocity. It is particu¬ 
larly challenging to solve such kinetic models, even us¬ 
ing advanced numerical techniques. In addition, a full 
kinetic description of the system might provide unnec¬ 
essary information, depending on the phenomena under 
investigation. This is the justification, e.g., for gyroki- 
netic theories where a strong magnetic field assumption^ 
leads to the perpendicular component of the velocity be¬ 
ing replaced by the magnetic moment considered as an 
adiabatic invariant. Similarly, anisotropy due to a strong 
magnetic field is also used to reduce the complexity of the 
original kinetic problem in the double adiabatic theory^. 
More generally, and ideally, reduced models obtained by 
some kind of closure leads to a decrease in the complexity 
of the original kinetic problem, while maintaining accu¬ 
racy and providing physical insights about the processes 
at work. Consequently, fluid reductions of kinetic equa¬ 
tions are often sought. 

Generally speaking, fluid models are obtained by pro¬ 
jecting the distribution function as follows: 

N 

/(x,v,i) ~ y^Pi(x,i)ei(v), (1) 

i=0 

where Pi(x, t) is the Tth fluid moment defined as the i-th 
moment of the distribution function with respect to the 
velocity v, i.e., 

Pi(x, t) = J y (g) V - • (g) v/ (x, V, t) dv, 

i times 


for all i G N and {ei(v)} is some specific fixed set of 
basis functions. The dynamics of the distribution func¬ 
tion /(x, V, t) is then given by the dynamics of the fluid 
moments that are functions of the configuration space 
coordinate x only. This makes fluid models, which in¬ 
volve quantities such as, the density p(x, t) = Po(x, t), 
the fluid velocity u(x, t) = Pi (x, t )/ Pq (x, t) and the pres¬ 
sure P(x, t) = P 2 (x, t) — Pi (x, t)/Po(x, t), convenient to 
interpret. Furthermore, since fluid variables only depend 
on X at each time, they are substantially less expensive 
to solve numerically than their kinetic counterpart. 

Clearly accurate reduced fluid models are desirable, 
but finding effective fluid closures is a difficult and largely 
open problem. Indeed, despite their strong physical rele¬ 
vance, a general or optimal procedure for obtaining them 
for the Vlasov equation does not exist. For example, con¬ 
sider the following simple free advection equation: 

5t/ = -v-V/, (2) 

which is the Vlasov equation with field dynamics re¬ 
moved. Multiplying Eq. ([2j) by and integrating with 
respect to the velocity, yields the following infinite hier¬ 
archy of moment equations: 

atPi = -v-p,+i, (3) 

for all i G N. In order to be able to solve Eqs. m , one 
has to truncate the infinite set of equations at some order 
V G N. However, because the time evolution of P/v in¬ 
volves Pat+i, the latter must be neglected or expressed in 
terms of lower order moments, i.e., Pat+i = Pat+i (Pi< at)- 
This is the ubiquitous closure problem for fluid reductions 
of kinetic equations. 

In conventional fluid closure theory a collision pro¬ 
cess and assumption of local thermodynamic equilibrium 
is basic. Instead, in this article we consider Hamilto¬ 
nian fluid reductions, where we investigate closures based 
on whether or not they preserve Hamiltonian structure. 
This allows us to select out a subset of all possible fluid 
closures that preserve the geometrical structure and this 
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prevents the introduction of non-physical dissipation in 
the resulting fluid moment system, without requiring 
nearness to thermal equilibrium. 

The case of two moments^ corresponds to the well- 
known exact water-bag reduction, so it is not surprising 
that it is Hamiltonian. However, as one increases the 
number of moments, this increases the dimension of the 
subset Pi<N^ and the constraints needed to preserve the 
Hamiltonian structure become more difficult to solve^i^, 
so that eventually, it is not possible to obtain a general 
analytic expression for the closure P/v+i = PN-\-i{Pi<N)‘ 

The problem of deriving Hamiltonian fluid models can, 
however, be tackled from different angles. Indeed, in¬ 
stead of Eq. other representations of the distribu¬ 
tion function can be used to decrease the complexity of 
the initial problem. In this paper we consider a gen¬ 
eral water-bag model^“— , which has also been used, e.g., 
in gyrokinetics^Tii^. In one dimension, this projection is 
obtained by replacing the distribution function with a 
piecewise constant function in the velocity v such that 

AT+l 

/(a;,?;,t) ~ y] ai©[v-Vi(a;,t)], (4) 

where are constants, 0 denotes the Heavyside distribu¬ 
tion and is a set of contour velocities. Like with 

the fluid moment projection, the dynamics of the dis¬ 
tribution function f{x^v^t) defined on phase space has 
been replaced by the dynamics of + 1 fields defined on 
configuration space, namely Vi(x, t) for all 1 < i < N-\-l. 

The use of the water-bag projection constitutes an ex¬ 
act reduction and consequently the resulting system is 
intrinsically Hamiltonia n^^d^ . When the number of field 
variables expressed in terms of fluid moments and the 
number of contour velocities are the same, the water-bag 
projection is easier to handle than that of the usual fluid 
moments representation; in particular, this is the case 
for the computation of a Poisson bracket. However, even 
though the contour velocities v^ are rather convenient to 
handle from a computational point of view, their macro¬ 
scopic physical interpretation is less obvious than for the 
fluid moments. Consequently, there is a balance to seek 
between the computational simplicity of the closure pro¬ 
vided by the water-bag model and the physical relevance 
of the fluid moments. 

In this article, we investigate links between the water- 
bag and the fluid moment representations in order to gen¬ 
erate new Hamiltonian closures. Indeed, any truncation 
of the infinite series given by Eq. (|lj) is preserved by the 
dynamics and hence constitutes a closure. As a conse¬ 
quence, the subset of all the water-bag distribution func¬ 
tions is invariant. Eollowing the water-bag projection, we 
perform a fluid reduction of the distribution function to 
obtain a Hamiltonian fluid model. Then, we construct a 
systematic procedure to obtain a fluid reduction from the 
water-bag distribution function by preserving the Hamil¬ 
tonian structure of the parent kinetic model. We extend 
this procedure to build general A/'-fleld Hamiltonian fluid 
models with N — 2 internal degrees of freedom. 


In Sec. El we provide the Hamiltonian structure of the 
Vlasov-Ampere equations which constitute the parent ki¬ 
netic model. The Casimir invariants of the associated 
bracket are provided. We introduce the water-bag dis¬ 
tribution function and give the associated Hamiltonian 
structure. Some properties of the system such as in¬ 
variants are discussed. In Sec. ED we establish a link 
between the water-bag and the fluid models. This is 
done by exhibiting a peculiar set of fluid variables that 
allows us to make explicit the fluid closure correspond¬ 
ing to the water-bag model. We use the density and 
the fluid velocity to account for the macroscopic energy 
of the system and we propose suitable variables to take 
into account internal degrees of freedom coming from mi¬ 
croscopic phenomena. The Hamiltonian structure of the 
resulting equations is provided and their Casimir invari¬ 
ants are discussed. We also address the thermodynamic 
implications of the new variables. Lastly in Sec. ElDl 
new models are proposed to extend the results obtained 
from the water-bag model to more general distribution 
functions. This allows us to construct general V-field 
fluid models that describe plasmas with N — 2 internal 
degrees of freedom. 

II. THE HAMILTONIAN STRUCTURE OF THE 
VLASOV-AMPERE EQUATIONS AND THE WATER-BAG 
MODEL 

We investigate the dynamics of a one-dimensional 
plasma made of electrons of unit mass and negative unit 
charge evolving in a background of static ions. This sim¬ 
plified system contains essential difficulties of more com¬ 
plete dynamics. We assume vanishing boundary con¬ 
ditions at infinity in velocity v and periodic boundary 
conditions in the spatial domain of unit length. The 
time evolution of the distribution function of the elec¬ 
trons f{x^v,t) and the electric field E{x^t) is described 
by the one-dimensional Vlasov-Ampere equations, 

dtf = -vdxf + Edyf, ( 5 ) 

dtE = -J, (6) 

where E = E — j E dx and J = j — J j dx are the fluc¬ 
tuating parts of the electric field E{x,t) and the current 
density j(x, t) = — f vf{x,v,t) dv^ respectively. 

The Vlasov-Ampere model possesses a Hamiltonian 
structure for the distribution function of the electrons 
f{x,v^t) and the electric field E{x,t). The Poisson 
bracket acting on functionals E[f,E] isi^“— 

{F,G} = J f[d,Ffd,Gf - d.Gfd^Ff 

+ Fed^Gf - Ged^Ff] Axdv, (7) 

where Ef and Ee denote the functional derivative of E 
with respect to the distribution function f{x,v,t) and 
the electric field E{x,t)^ respectively. Bracket ([7j) is a 
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Poisson bracket, i.e., it satisfies four essential properties: 
it is linear in both its arguments; it is alternating, i.e., 
{F,F} = 0; it satisfies the Leibniz rule, i.e., {F,GH} = 
{F, G}H + G{F, iL}; it verifies the Jacobi identity, i.e., 

{F, {G, H}} + {H, {F, G}} + {G, {H, F}} = 0. 

for all functionals F, G and H. The Hamiltonian of the 
system, which corresponds to the total energy, is 

= J y/ dxdv + j ^ dx, (8) 

where the first term accounts for the kinetic energy of the 
electrons and the second term corresponds to the energy 
of the electric field. Bracket ([7j) and Hamiltonian ([8|) lead 
to Eqs. and dni) by using dtf = {/, F} and dtE = 

{E^n}. 

The Casimir invariants of a bracket {•, •} are particular 
observables G that commute with all observable F, i.e., 
{F, G} = 0 for all functionals F. Bracket ([7j) possesses a 
local (x-dependent) Casimir invariant given by 

Cioc = d,E + J f dv, (9) 

which corresponds to Gauss’s law. There are also global 
(x-independent) invariants. Namely, 

E = J E dx, (10) 

which expresses the fact that the mean value of the elec¬ 
tric field remains constant. This results from the periodic 
boundary conditions in space and from the definition of 
the electric field F = —dx^ where ^(x,t) is the elec¬ 
trostatic potential. Finally, there is a family of global 
invariants given by 

Ci= J (pif) dxdv, ( 11 ) 

where 0(/) is any function of /. This family of Casimir 
invariants arises from particle relabeling symmetry and 
includes, e.g., the cases of conservation of the total mass 
and the usual entropy. 

The water-bag model is a particular solution of Eqs. 
and m with a piecewise constant initial condition for the 
distribution function /(x,'C,t), 

AT+l 

fN{x,v,t) = '^aiQ[v-Vi{x,t)], (12) 

i=l 

which can be done for any N E N. In water-bag theory 
one is interested in approximating a smooth initial con¬ 
dition by a water-bag approximation, such as that shown 
in Fig. [H In order for this distribution function to have 
compact support, we further require the following con¬ 
straint: 

AT+l 

y; fli = 0. 

i=l 



Figure 1. Sketch of a distribution function (in light red) and 
its water-bag approximation (in dark blue). 


Moreover, the velocities are supposed to be or¬ 

dered such that for all (x,t) G [0;l[xM+ we have 
vi(x,t) < V 2 (x,t) < ••• < VAr+i(^,t). In what follows, 
we will refer to this distribution function as an A/"-water- 
bag distribution function. A distribution function of the 
form of Eq. m is solution of Eqs. © and and hence 
its form is preserved by the dynamics, if and only if the 
contour velocities satisfy^”— 

dt^i = -^idx^i - E, ( 13 ) 

for all 1 < i + 1 and 

1 ^+1 

dtE = ( 14 ) 

i=l 

As a consequence, solving the Vlasov equation for a 
water-bag distribution function fjsj{x^v^t) is equivalent 
to solving the N -\-l contour equations given by Eq. m 
for the contour velocities Vi(x,t). Coupling between the 
different contours is then provided by Eq. ([TT|^ . 

The water-bag model possesses a Hamiltonian 
structurei^i^Si^ inherited from the original Vlasov- 
Ampere equations. This means that there exist a bracket 
{Awb and a Hamiltonian F such that Eqs. ([T^ - (p!T|) are 
obtained by dt^i = for all 1 < i < A/" H- 1 and 

dfE = {F, respectively. For the water-bag dis¬ 

tribution function of Eq. m, the dynamical variables 
are the contour velocities Vi(x,t) and the electric field 
F(x, t). Using the chain rule, the functional derivative of 
F with respect to v^, denoted Fi, is 

F, = ( 15 ) 

for all 1 < i < V + 1. Along with Eqs. ([12]) and Eq. (fT^ 
this leads to the water-bag bracke t 


{F, G}wb 



-FAGi + GiFs 

di 


FiGs 


dx. 

( 16 ) 
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One can show that Bracket m is a Poisson bracket, 
which is a property inherited from the Vlasov-Ampere 
equations. The Hamiltonian of the system is 


III. LINK BETWEEN THE WATER-BAG MODEL AND 
THE FLUID MOMENTS OF THE DISTRIBUTION 
FUNCTION 





which is obtained from Eq. (|H)) by using Eq. (HU). 

An important feature of Bracket m is that it is 
closed. Thus, for any number of bags V G N, the set 
of the water-bag distribution functions /at of N bags is a 
sub-Poisson algebra of the Vlasov-Ampere model. This 
means that the Vlasov-Ampere dynamics preserves the 
number of bags. In particular, the water-bag model is 
Hamiltonian for any number of bags. This is particularly 
interesting from a numerical point of view, e.g., as there 
is no nonphysical dissipation introduced by the water-bag 
approximation even for a low order approximation with 
a small number of bags. 

Bracket ([T6|) possesses several Casimir invariants. By 
using Eq. ([T2j), Eq. (|9]) becomes 


AT+l 

^loc ~ 9xE ^ ^ 
i=l 

The global invariant given by Eq. m is preserved by 
Eq. ([7j)- The family of Casimir invariants given by 
Eq. (fTT)) is projected to 


N 

Cl ='^(j){Ai) (vi+i - Vi), 

i=l 


The contour velocities Yi{x,t) provide immediate ki¬ 
netic theory information: they define the partitioning 
of the distribution function in the velocity space, sort¬ 
ing particles into water-bags according to their velocities. 
However, their interpretation on the fluid level in terms 
of moments Pi(x,t), given by 

Pi(x,t) = JvY{x,v,t) du (18) 

for all i G N, is not so clear. This relationship is given 
explicitly by inserting Eq. (|12p into Eq. m, yielding 

_ AT+l 

Pi{x,t) = akvl+^{x,t), (19) 

^ ^ k=l 

for all i G N. Our strategy is to use the information 
provided by the water-bag model to build Hamiltonian 
models for the fluid moments P^(x,t). Indeed, as stated 
in Sec. im the water-bag model possesses a Hamiltonian 
structure. As a consequence, by expressing the contour 
velocities with respect to the fluid moments we can ob¬ 
tain particular fluid models with an arbitrary number 
of moments. We believe that this is a useful strategy 
because constructing general fluid models can be tech¬ 
nically very challenging^. The closures provided by the 
water-bag models provide insight for building more gen¬ 
eral fluid models. 


A. The single water-bag model 


where Ai = and 

k=l 



for all 1 < i < A" + 1. As Cl is a Casimir invariant 
for any function 0, this shows that the projection of the 
invariant given by Eq. dl) leads to the generation of N 
invariants, namely 

Ci^f = v^-j-i v^, 

for all 1 < 2 < A. However, Ci is computed such that 
[f,6Ci/6f] = 0 for any distribution function /, where 
[g^h] = dxgdyh — dxhdyg. If we now look only at water- 
bag distribution functions, the requirement for Ci be¬ 
comes [f,6Ci/6f] = 0 for all / given by Eq. ([I2]), and 
hence is less restrictive. This leads to the creation of 
an additional invariant, e.g., vi. Thus, there are N 1 
Casimir invariants given by v^ for all 1 < i < A -h 1, i.e., 
as many Casimir invariants as the number of fields ^2^ 


Eirst consider the case of a single water-bag in order to 
illustrate our approach, which will be generalized to an 
arbitrary number of water-bags corresponding to an arbi¬ 
trary number of fluid moments. This simple model con¬ 
stitutes a good illustration of our strategy. If we consider 
a single water-bag, or equivalently two contour velocities, 
the distribution function simply reads 

fi{x,v,t) = 0[i; - Vi{x,t)] - ©[v - V 2 {x,t)]. 

where we have set ai = 1 defining unit height to the 
water-bag. Using Eq. m, the first two moments of the 
distribution function are 

Vo — V? 

Po = V 2 - vi and Pi =- - -. 

which upon inversion yield 



Defining the density p = Pq and the fluid velocity u = 
Pi/Po, we obtain the fluid variables 

Pa , ^ 

Yi = u - - and Y 2 = uP - . 
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As a consequence, we are able to express the contour 
velocities with respect to the usual fluid moments. In 
terms of these fluid variables, the water-bag bracket given 
by Eq. is 

{F, G}i = J [Gud.Fp - Fud.Gp + GuFe - dx, 

where now Fp and Fu denote the functional derivative 
of F with respect to p and respectively. This bracket, 
which corresponds to the cold-plasma bracketii, is closed, 
a property inherited from the original water-bag bracket 
given by Eq. dini). In terms of the variables p, u and E, 
Hamiltonian (fT7|) becomes 

'H[p,u,E] = j da;, 

where U{p) = p^/24 is the specific internal energy of the 
system. The pressure is defined by the usual thermody¬ 
namic relation P{p) = p^dU/dp = p^ll 2 . The reduced 
moments, defined by 

Siix,t) = J{v- uyf{x,v,t) dv, (20) 

for alH > 2, appear to be suitable variables to describe 
the Poisson structure of the fluid equations resulting from 
the Vlasov-Ampere model^. With this definition, the 
second order reduced moment reads S 2 = P/p^^ which 
eventually leads to S 2 = 1/12, i.e., S 2 is constant. (Note, 
throughout we will express the Si in l/a{ units for all 
i > 2.) 

In term of the reduced moments, the Hamiltonian of 
the system reads 

n[p,u,E] = ^ J {pu^ + p^S 2 + E^) dx, 

and the equations of motion are 

dtp = {p,'H}i = -dx{pu), (21) 

dtu = {u, H}i = -ud^u - -dx{p^S 2 ) - E, (22) 

P 

dtE = {E,n}i = p^. (23) 

These are the equations for a barotropic fluid undergo¬ 
ing an adiabatic process. Indeed, the relationship be¬ 
tween the pressure P and the density p is such that 
P/p^ = S 2 = 1/12 is a constant. This is the charac¬ 
teristic isentropic equation of state of an ideal gas with 
one degree of freedom. The local Casimir invariant given 
by Eq. dH) is preserved and is 

C\oc = dxE + p. (24) 

In addition, the system have three global invariants: the 
mean value of the electric field given by Eq. m and 


p = J pdx, 

(25) 

u = u dx. 

(26) 


The first invariant is conservation of total mass, which 
results from the fact that the system is isolated. This 
Casimir invariant results from the projection of Eq. ([TT]) . 
As for the electric field, the last Casimir invariant, corre¬ 
sponds to conservation of the mean value of the velocity. 
As noted, this additional invariant arises from the closure 
procedure. 


B. Two water-bag model: introduction of the 
thermodynamical variables 


Now consider the case of two water-bags, in order to 
characterize more precisely the closure provided by the 
water-bag model and its relation to the fluid moments. 
Indeed, there exists an infinite number of Hamiltonian 
fluid models with three moments^. By using the reduced 
moments defined by Eq. (EUD . Hamiltonian models for 
the variables p, u and S 2 are such that S 3 is an arbitrary 
function of 5'2- Two water-bags means we have three 
contour velocities or, equivalently, three fluid moments, 
which provides a particular example of the more general 
closure S 3 = 5'3(5'2). Such a distribution function, whose 
expression is given by Eq. da with V = 2, is represented 
in Eig. [21 By using Eq. ([T9|) . the first three moments of 
the distribution function are 


Pb = V3 - vi + a2(v3 - V2), 


Pi = 


+ CL 2 ■ 


P 2 = 


Ci2~ 


Vq — Vo 


Unlike the single water-bag case, expressing the contour 
velocities with respect to the fluid moments with two 
water-bags is more complicated. This is partially due 
to the fact that P 2 involves cubic terms in the contour 
velocities v^. This issue becomes more acute as the num¬ 
ber of water-bags increases. A first step toward the fluid 
representation can be easily done by using the following 
variable: 


_ _ V2 - Vl _ 

V3 - Vi -h a2(v3 - V2) ’ 

along with the density p, the fluid velocity and the 
electric field E as used in the single water-bag model. In 
this case, pni simply corresponds to the density of the 
particles contained in the first bag of unit height. We see 
that J pni dx, i.e., the total number of particles in the 
first bag, is a Casimir invariant of Eq. ([T6|) . hence is a con¬ 
stant of motion. The contour velocities can be expressed 
explicitly with respect to these variables as follows: 


Vl = + - 


Y 2 = UP 


ys = Up 7: 


pa2ni{ni - 2 ) - 1 

2 l + a2 ’ 

pni{a2ni + 2 ) - 1 

2 l + a2 ’ 

pa2nl -h 1 
2 1 -j- a2 
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Expressed in terms of the variables (p, ni, E), ([T^ 
takes the particularly simple form 


{F,G}2 = I 


GudxFp — FudxGp + GuFe — FuGe 


- -{F,Gu - GiFja.m + 

P Ci 2 p 



dx, 

(27) 


where Fi denotes the functional derivative of F with re¬ 
spect to Til. Bracket (|?7|) is closed which is, as in the 
single water-bag model, a property inherited from the 
original water-bag bracket given by Eq. (fT6j) . In terms of 
the variables p, rii and E, Hamiltonian (HH becomes 


'H[p,u,ni,E] = ^ j [pv? + S 2 {nx) + dx, (28) 
where 


S2{ni) 


1 6 a 2 n^ + 4a2(a2 — l)n\ — Sa^nf 

12(l-ha2)2 


(29) 


This shows that *52 is a function of rii only. The spe¬ 
cific internal energy now becomes U (p, ni) = p^5'2(^i)/2. 
As a consequence, in addition of the pressure defined 
as the thermodynamic conjugate variable of the density 
through the relation P = p^dU/dp^ we can define some 
potential pi as the conjugate variable of ni, such that 
Pi = p^S 2 {ni)/2. This shows that more accurate fluid 
models are obtained by introducing more information on 
the thermodynamic properties of the system through in¬ 
ternal degrees of freedom. Bracket ([27j) and Hamilto¬ 
nian ([28]) lead to Eqs. ([2T])-([23|) and the following addi¬ 
tional equation: 


dtm = {ni,'H)2 


—udxTii + 


1 + a2 
a 2 p 


dx 



The first term of this equation is an advection term, while 
the second is a flow term resulting from the potential 
Pi = P^5'2(ni)/2. Indeed, analogously to the pressure P 
that drives a force —dxP in Eq. (E21), here pi drives a 
flow dxPi- 

Along with the Casimir invariants given by Eqs. 
m, and the two water-bag model has the following 
global invariants: 

pni dx, (30) 

The Casimir invariant given by Eq. (|5n|) is inherited 
from the original Vlasov-Ampere model and amounts to 
conservation of the total entropy. Indeed, in statistical 
physics, the entropy of a system is related to its num¬ 
ber of microstates. Here the microstates are given by the 
amount of particles in each bags. The invariant given by 



Eq. ([31]) is a new conserved quantity and is generated by 
the reduction procedure. Noting that Eq. ([29|) defines a 
bijection g : rii e [0; 1 ] S 2 = p(^i) G M+ such that 
we can write rii = ^{ 82 ) where k = the previous 
invariants become 


Pk{S 2 ) = J pf^{S 2 ) dx. 

Therefore, we see that the two water-bag model is a 
particular case of the more general three moments fluid 
model^. 

As stated previously, the two water-bag model is 
closed. Thus, we expect the associated fluid model to 
be closed too, and the fourth reduced moment S3 to be 
a function of S 2 only. In what follows, we choose a 2 > 0. 
This corresponds to a configuration in which the second 
bag is taller than the first, as depicted in Eig.[ 2 l The case 
a 2 < 0 is equivalent through the symmetry v ^ —v and 
rii 1 — rii. By using the definition of the third order 
reduced moment given by Eq. (|2Q|) with i = 3, we find 

e _ a2{ni - + a2nif 

4(l + a2)3 

Since S3 is a function of ni only and rii = 1 ^( 82 ), we see 
that S3 is a function of 82 only, i.e., S3 = 83 ( 82 ) as ex¬ 
pected. This relation was expected because it was shown 
in Ref. 0 that it is the case for general Hamiltonian clo¬ 
sures with three fluid moments obtained from the Vlasov 
equation. The fraction of particles in the first water-bag 
parametrizes the curve ni ^ [S 2 (ni), 83 ( 711 )]. This re¬ 
sult corresponds to a closure for the heat flux q = p^ 5 ' 3/2 
as a function of the density p and the pressure P = p^ 82 - 
We do not give the explicit relationship between 82 and 
S3 here because it does not provide much information. 
However, the dependence of S3 on 82 is plotted in Eig.[ 2 ] 
for different values of a 2 . We observe that as a 2 —> 0, we 
have 82 1/12 and > 0. This is consistent with the 

results of Sec. lHI Al Indeed, for a 2 = 0 the two water-bags 
have the same heights and, as a consequence, merge into 
one such that we recover the values of the fluid moments 
corresponding to a single water-bag model. Moreover, as 
a 2 increases, the solution tends rapidly toward an equi¬ 
librium obtained with a 2 + 00 . 

In summary, the two water-bag model can be conve¬ 
niently described by using appropriate fluid variables: 
the density p and the fluid velocity u are natural fluid 
quantities that take into account, in particular through 
the definition of the kinetic energy K = f pvS j 2 dx, the 
macroscopic energy of the system. The internal effects 
are described by the internal energy. Using the parti¬ 
tioning of the number of particles into two bags as the 
internal degree of freedom of the system, we can define 
the specific internal energy as I/(p, ni) = p^ 5 ' 2 (^i)/ 2 . 
The two water-bag model corresponds to a system with 
one internal degree of freedom described by ni. 
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the particularly simple form 


{F,G}3 = j 


GudxFp — FudxGp + GuFe — F^Ge 


-{FiGu - GiF^)d^ni + ( — 


dx, (32) 


where Fi denotes the functional derivative with respect to 
rii for i G {1,2} and where the summation over repeated 
indices from 1 to 2 is assumed. Here yd is a constant 2x2 
symmetric matrix given by 

^ (1 + Q 2 ) f as -as A 

a2as \~^3 F a2){a2as) J 


In terms of the variables p, ni, 77 - 2 , and Hamilto¬ 
nian (HU) becomes 


Figure 2. Upper panel: sketch of a double water-bag distri¬ 
bution function. Lower panel: plot of S 3 as a function of S 2 
for a double water-bag distribution function corresponding to 
the upper panel and for different values of a 2 (given by the 
colorbar). 


C. Three water-bag model 

In this section, we demonstrate the usefulness of the 
new thermodynamical variables for linking the water-bag 
and fluid models by considering the three water-bag dis¬ 
tribution function, whose expression is given by Eq. m 
for N = 3. A three water-bag model is equivalent to a 
Hamiltonian fluid model with four fluid moments. Even 
though a particular closure based on dimensional anal¬ 
ysis has been previously found for such fluid models^, 
there is currently no general Hamiltonian closure for fluid 
models with four moments. Einding all the closures 
is difficult because it requires solving non-linear PDEs 
obtained from the Jacobi identity. We show here, by 
using the thermodynamic variables, that the water-bag 
distribution function with three water-bags is another, 
parameter-dependent closure for fluid models with four 
moments. This solution is simpler to compute than the 
general closure, and gives us some useful information. 
Different three water-bag distribution functions are de¬ 
picted in Eigs.[3l[l]and[5l 

Eollowing the procedure of Sec. MbI we introduce the 
thermodynamic variable 


'H[p,u,ni,n2,E] = ^ J [pu^ + 82(711, n2) + E'^] da;, 

(33) 

where 5'2(^i,^2) is given by Eq. (|37|) . As in the two 
water-bag model, we can define the thermodynamic po¬ 
tential fj, = dU/dn = p‘^{dS 2 /dn)/ 2 . Analogously to 
multi-components systems, the variables rii act as in¬ 
ternal degrees of freedom that characterize the model 
through a partitioning of the particles with respect to 
their energy or, equivalently, their temperature. In¬ 
deed, it is expected that a system like the collisionless 
Vlasov-Ampere system should not thermalize and as a 
consequence is described by more than one temperature. 
Bracket (Ea) and Hamiltonian El lead to Eqs. EH - 
El, and in addition the following equations: 

— {^2,^}3 — T l^ik^x 

for i G {1,2}, where summation over repeated indices 
form I to 2 is assumed. These equations exhibit an ad- 
vection term and a driving term through the existence of 
a potential /j, = p‘^{dS 2 /dn)/ 2 . 

Along with the Casimir invariants given by Eqs. El, 
d, El, and El, Bracket (|3^ has the following global 
invariants: 


V 2 drik) ’ 


pn 2 


pn2 da:, 


G 2 



a2n‘l 
(1 +a 2 ) 


asjTii + n2)^ 

(1 + a 2 )(l + 02 + 03 ) 


} 


dx . 


_ _ (1 + 02 )(V3 - V 2 ) _ 

V4 - Vi + a2{v4 - V2) + a3(v4 - V3) ’ 

in addition to p, tx, ni, and E used in the two water- 
bag model. Here pn 2 corresponds to the density of the 
particles contained in the second bag. The expressions of 
the contour velocities in terms of these variables are given 
by Eq. (|M|) and will not be detailed here. Expressed in 
terms of the variables (p, xx, rxi, 7 x 2 , E), Bracket ([TB|) takes 


Thus, as for the other models, we have as many global 
Casimir invariants as dynamical variables. Moreover, 
there is some generalized velocity which is a common 
feature of all fiuid models derived from Vlasov-Ampere 
equations. 

By using Eq. EZD, we can express £' 2 , Ss and S 4 , 
as functions of rxi and 7 x 2 - Eor any (a 2 ,a 3 ), these 
functions define a unique two-dimensional manifold 
[S' 2 (ni, n 2 ), S' 3 (ni, n 2 ), 8 ^( 711 , 712 )]. 
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Figure 3. Upper panel: sketch of a three water-bags distri¬ 
bution function whose typology exhibits a bell-shape. Lower 
panel: colormap of 5^4 as a function of aS '2 and S 3 for (a 2 , ns) = 
(2,—1.75) corresponding to a distribution function given by 
the upper panel. 


There are mainly three configurations of interest for 
the distribution function, with other typologies obtained 
by using different symmetries. The first configuration has 
0 < —as < a 2 . This bell-shaped configuration is shown 
in Fig. [3l The second configuration has 0 < —a 2 < as 
and corresponds to a case in which the second bag is 
the smallest. Such a typology exhibits a “hole" in the 
distribution function shown in Fig. 01 The last typology 
has a 2 > 0 and as > 0 and corresponds to a configuration 
with the third bag taller than the second, which is taller 
than the first. Such a typology is shown in Fig. [5l In 
Figs. [3l m and [5] we plot *54 as a function of S 2 and S 3 
for the three configurations. 

Observe, despite the change in distribution function 
typology, which result in a change in the typology of the 
manifolds defined by the closure, *54 always increases as 
S 2 increases. Moreover, for a configuration as depicted 
in Fig. [5l the sign of S 3 is fixed, whereas it may vary 
for the other configurations depending on the respective 
widths of the water-bags. The same computation can 
be performed for higher order moments. In particular, a 
Hamiltonian fluid model for four moments requires clo¬ 
sures on the fourth and fifth order moments respectively, 
namely S 4 and 83 ^. 


D. N water-bag model 

The method presented in the previous subsections can 
be extended to an arbitrary number of water-bags with a 



0.04 0.06 0.08 0.1 


^2 

Figure 4. Upper panel: sketch of a three water-bags distri¬ 
bution function whose typology exhibits a hole. Lower panel: 
colormap of aS '4 as a function of aS '2 and S3 for (a 2 ,a 3 ) = 
(—0.15, 0.9) corresponding to a distribution function given by 
the upper panel. 

corresponding arbitrary number of fluid moments. This 
is important because by increasing the number of water- 
bags we increase the accuracy of the description, allow¬ 
ing for more refined kinetic effects. This is consistent 
with the fact that the water-bag models come from a dis¬ 
cretization of the distribution function in velocity space. 
An analogy can be made with vibrations of structures^. 
In these systems, the frequency spectrum is continuous. 
However, these models can be accurately described by a 
finite number of coupled springs with a discrete spectrum 
as long as the number of springs is sufficiently high. Thus 
by increasing the number of water-bags, yet keeping it fi¬ 
nite, we can recover kinetic information about the system 
within the framework of a fluid description. Considering 
N bags, we define the following variables: 

AT+l 

P=-^ 

i=l 

AT+l 

E 

i=l 

^ ~ AT+l ’ 

2 

k=l 

I 

(v^+l - Yl) 

~ AT+l ’ 

Z] a/eV/c 
k=l 

for all 1 < / < A" — 1. In this case, prii corresponds to 
the density of the particles contained in the i-th bag for 
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x10'^ xIO'^ 



Figure 5. Upper panel: sketch of a three water-bags distribu¬ 
tion function whose typology exhibits a monotonic increase. 
Lower panel: colormap of S 4 as a function of aS '2 and S3 for 
(a 2 ,a 3 ) = (0.25,2) corresponding to a distribution function 
given by the upper panel. 


1 < i < N — 1. We see that f prii dx, i.e., the num¬ 
ber of particles in the Uth bag, is a Casimir invariant of 
Eq. (pU|) . The contour velocities can be expressed explic¬ 
itly with respect to these variables such that 

Y- =u^ p^i (34) 

for all 1 < i < A/" + 1 and where 

N-l N -1 1 2 

1 - 2^ nfe - fflAT+i 2^ — , 

k=l k=i 


N 




AT+l 




= 1 ^AT+l 


and 


/3 is di constant tridiagonal (A^ — 1) x (A^ — 1) symmetric 
matrix given by 

/ Ai -Ai 0 ... 0 \ 

—Ai Ai + A2 ~A2 ... 0 

p — 0 —A2 A2 + A3 ... 0 

Vo 0 0 ... Xn-J 

where 

i i-\-l 

, _ k=l 1=1 

Ai — 

Bracket (|35]) can be further simplified. Indeed, noting 
that p is symmetric, hence diagonalizable, we introduce 
the variables 


Vi = 

fe=l 

for 1 < i < N — 1. The quantity f pVi dx corresponds 
to the cumulative number of particles in the i first bags. 
Eventually, Eq. (|35]) takes the even simpler form 


{E,G}w = 


1 


/ 


GudxFp — FudxGp + GuFe — FuGe 


--{FiGu - GiF^)d,Vi + XiX±d, ( U 


dx, (36) 


where Fi denotes the functional derivative of F with re¬ 
spect to for 1 < i < A/" — 1 and where summation 
over repeated indices from 1 to A^ — 1 is assumed. The 
Hamiltonian associated with this model is given by 


= \ j + P^S 2 {vi,...,vn-i)+E‘^] da;, 

where the reduced moments can be computed from 
Eq. ([2Q|) and are given by 


= ^AT+l + 


1 


aw+i 


N-l 


k=l 


N-l 

ajv +1 ^ 

k=m 


rik 

Ak 


Si{x,t) 


(j + l)pi+l E 


for all 1 < m < A/". Expressed in terms of the variables 
(p, iz, ni,..., nyr-i, E), Bracket ([16]) takes the particu¬ 
larly simple form given by 


{F,G}n = 


/ 


GudxFp — FudxGp + GuFe — FuGe 


-{FiGu - GiFu)d,ni + pik-d, (^ 
P P \ P 


dx, (35) 


where Fi denotes the functional derivative with respect 
to for 1 < i < A/" — 1 and where the summation over 
repeated indices from 1 to A^ — 1 is again assumed. Here 


for all i >2. By using Eq. (|34|) . this eventually becomes 

AT+l 


where 


= 


N 


i V “P J. 

(j + 1) E akCk^'^ivi,...,VN-i), (37) 


?JV +1 = \^ 

2 = 1 

Cm ~ Cat+I T 


ai 


AT-l 


= 1 ^AT+l 

1 


aw+i 


~ ^k-l 

i — VN-I — CLN +1 / ^ 
k=i 
AT-l 

1 — Z/AT-l — aw+l ^ 


^k 

k'k — k'k-l 


k=m 


Ak 
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Figure 6. Closures of the two water-bags model that define 
the projection of the edges of the closure of the three water- 
bags model. 


^11 lilfTf Tl llTrtT 


V 



for all m such that 2 < m < N — 1. This shows that the 
reduced moments Si are functions of the thermodynamic 
variables Ui only. The equations of motion of the system 
are given by Eqs. and 


/ ovx fp^dS2\ 

otVi = {vi,n.]N = -udxVi + —Ox I ) > 

for all 1 < z < — 1. Along with the Casimir invariants 

given by Eqs. (IM . (fTUl) . and (123), the N water-bag model 
has the following global invariants: 



dk 


a 


2 

AT+l 




AT-l 

+ Uat+I 

l=k 



dx. 


for all 1 < i < — 1. 

The closure provided by the N water-bag model is not 
straightforward. Indeed, the closure is such that S'at+i = 
5'Ar+i(S'2,..., S'at), which defines an — 1-dimensional 
manifold in parametrized by (z/i,..., z/tv-i)- How¬ 
ever, one can apply the tools developed throughout the 
previous sections to visualize such a manifold. Indeed, 
consider, e.g., a three water-bag distribution function as 
shown in Fig. [5] and consider all the distribution func¬ 
tions (three in this example) obtained by combining two 
of these bags. These are shown in Fig. [6l We see that 
the edges of the manifold defined by the four fluid mo¬ 
ment closure S 4 = 84 ( 32 ^ Ss) correspond to the closures 
of the three fluid moments models associated with every 
combination of two water-bags of the initial three water- 
bags distribution function. Analogously, the projections 


Figure 7. Upper panel: sketch of a bell-shaped water-bag 
distribution function with twenty-seven bags. Lower panel: 
projection of the edges of the manifold dehning the closure 
for the distribution function given by the upper panel. 


on the (5'2,5^3, *54) space of the edges of the manifold 
defined by the closure Sn-\-i = *5Ar+i(*52,..., aSat) in the 
(*^ 2 , S's,..., ^AT+i) space correspond to the closure of the 
three fluid moment models associated with every com¬ 
bination of two water-bags of the initial N water-bag 
distribution function. 

We illustrate the above edge description with the fol¬ 
lowing example. Consider a Maxwellian distribution ap¬ 
proximated with a twenty-seven water-bag distribution 
function as shown in Fig. [71 The corresponding fluid 
closure is such that S 28 = *^ 28 (* 52 , • • •, *527)^- The projec¬ 
tion of this high dimensional manifold on the (*52, *53, * 54 ) 
space is depicted in Fig. [71 Consequently, we see that 
the information about the whole system is given by all 
the possible couplings between two different water-bags. 
This makes the study of systems with a high number of 
fields easier as it eventually reduces to the study of cou¬ 
pled subsystems with three fields. 

Inserting arbitrary functions in front of the terms 
{Fi/p)dx{Gi/p) in ([36|) may allow us to extend this 
bracket to more general Poisson brackets of the form 


/{ 


{F, G}= \ GudxFp - FudxGp + G^Fe 


— FuGe - {FiGu — GiFu)dxk'i 

P 


+ cri(z/i) 



(38) 


where the are arbitrary functions. In addition, we 
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consider Hamiltonians of the general form 

'^ = \j [pu^ + P^S‘ 2 ,{vi,...,vn-i)+E‘^] da;, 

where now the dependence of S 2 on is arbitrary. The 
choice of the unknown functions (Ji and the reduced mo¬ 
ment • • •, ^AT-i) should be based on physical ar¬ 

guments. As an example, we expect homogeneous ini¬ 
tial conditions (i.e., p = po, '^ = 0 and Vi = z/^o) to 
be linearly stable. This leads to constraints on cji and 
• • •, It is possible show, e.g., that the two 

water-bag closure is always stable with respect to a ho¬ 
mogeneous equilibrium. Along with the ones given by 
Eqs. ([25]) and ( 1101 ), the global Casimir invariants of the 
extended water-bag bracket given by Eq. dSSl) are 



where z^J = l/y^ for all 1 < I < N — We also 
notice that Brackets (|38|) for A" = 2 are the most general 
Poisson brackets with three moments^. Whether or not 
this is the case for any N is an open question. 


IV. SUMMARY 
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i=l 


In summary, we exhibited a method for constructing 
Hamiltonian fluid models with an arbitrary number of 
fluid moments from the Vlasov-Ampere system. This 
construction relies on the Hamiltonian structure of the 
water-bag representation of a distribution function. We 
introduced suitable fluid variables, based on thermody¬ 
namic considerations, to replace the less meaningful con¬ 
tour velocities. The density and the fluid velocity were 
used to describe macroscopic phenomena, while the par¬ 
titioning of the particles into the different bags was used 
to define internal degrees of freedom in the system, ac¬ 
counting for microscopic effects. By using these variables, 
we were able to link the water-bag and fluid models and 
to make explicit the corresponding closures. We showed 
that, for an arbitrary number of water-bags, the general 
associated closure can be constructed from knowledge of 
the couplings between all the other water-bags. Based on 
these results, we proposed a general N field fluid model to 
describe plasmas with N — 2 internal degrees of freedom. 


where the heights bi are constants. As a consequence, there are 
2N fields (resp. N parameters), namely v^{x,t) and v^{x,t) 
(resp. bi) for all 1 < -i < A. By using our notation, this for¬ 
mulation can be written fj^ix^v^t) = f 2 N-i{x^v^t) where ai = 
-a 2 N+i-i = k, yi{x,t) = vr{x,t) and V 2 N+i-i{x:t) = v+ {x,t) 
for all 1 < i < A. Gonsequently, the formulation given by this 
representation only has an even number of contour velocities 
Vi{x,t) (field variables) and only has A independent free pa¬ 
rameters. This corresponds to A adiabatic fiuids. 
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Indeed, a water-bag distribution function leads to the constraint 
AT+l 

S{v - Vi)dx^^ = 0 . 
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